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^ ' The possibihty of CP violation is studied in the Higgs sector of a supersymmetric standard 

^5 ■ the CP symmetry. We find that CP violation is viable in this model at the one-loop level, in an 



explicit way, if the radiative corrections from the third generation of quarks and squarks are taken 



,^ ' into account. In the presence of explicit CP violation, at the one-loop level, the upper bound on 

ji« ■ the mass of the lightest neutral Higgs boson and the productions of the neutral Higgs bosons via 

(D ; , _ 

r~| , the Higgsstrahlung process in high-energy e"^e collisions are calculated. We find that the upper 

bound on the mass of the lightest neutral Higgs boson increases as the number of Higgs singlets 

;_, ■ increases in a regulated manner. The production cross sections of the neutral Higgs bosons also 

show a reasonable increasing behavior with respect to the number of Higgs singlets. 
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I. INTRODUCTION 

The Standard Model (SM) possesses a complex phase in the Cabbibo Kobayashi Maskawa 
(CKM) matrix, but the size of the complex phase in the SM, as determined by several 
precision experiment, is too small to account, in terms of the Sakharov mechanism [1], for 
the baryon asymmetry of the universe [2]. In this context, a study has been devoted to 
the possibility of explaining CP violation in the minimal supersymmetric standard model 
(MSSM) [3]. The supersymmetry (SUSY) itself is widely regarded as an essential ingredient 
for the new physics beyond the SM [4] . Supersymmetrically-extended standard models may 
have various sources of complex phases, either explicitly or spontaneously. These complex 
phases may arise from the Higgs sector during the spontaneous symmetry breaking or from 
radiative corrections at higher level. Also, the soft SUSY-breaking terms, which are required 
as essential ingredients to construct the full Lagrangian density for practical analysis of the 
supersymmetric models, might be complex. 

The minimal supersymmetric standard model, which has no Higgs singlet and, thus, is 
the simplest version of supersymmetric standard model, has been investigated exhaustively 
with respect to various aspects of phenomenology [5], and the possibility of the CP violation 
[6]. It has been noticed that it can accommodate CP violation at the one-loop level, as 
complex phases in the soft SUSY-breaking terms may induce explicit CP mixings between 
the scalar and the pseudoscalar Higgs bosons in the neutral Higgs sector [7]. 

The Higgs sector of the MSSM may be enlarged by adding a Higgs singlet. This model is 
known as the next-to-MSSM, and might be referred to as the (M-l-l)SSM. One of the mo- 
tivations for the enlargement is to solve dynamically the so-called dimensional yU-parameter 
problem in the MSSM by means of the vacuum expectation value (VEV) of the Higgs sin- 
glet. The (M+1)SSM, including the CP phenomenology, has also been investigated by many 
authors [8-10]. At the tree level, the (M-l-l)SSM may contain at most one nontrivial CP 
phase after redefining the Higgs fields. At the one-loop level, it is observed that there are no- 
ticeable effects of explicit CP violation in the (M-l-l)SSM on the neutral and charged Higgs 
boson masses, when radiative corrections due to relevant particles and their superpartners 
are taken into account [10]. 

A further extension of the MSSM with two additional Higgs singlets is also possible. This 
model might be referred to as the (M-|-2)SSM. Actually, in the literature, some attempts 



have been reported to study the viabihty of spontaneous CP violation in the Higgs sector of 
the (M+2)SSM [11]. It has been found that, at the tree level, the spontaneous CP violation 
might occur in the Higgs sector of the (M+2)SSM. 

In principle, at the weak scale, the supersymmetrically-extended standard model, which 
has at least two Higgs doublets, may embrace additional Higgs multiplets. In other words, 
from the theoretical point of view, one may add an arbitrary number of Higgs singlets to 
the MSSM. Let us denote the MSSM with n additional Higgs singlets as the (M+n)SSM. 
In order for the (M+n)SSM to be theoretically interesting and consistent with experimen- 
tal constraints, its phenomenological implications should be studied. With respect to the 
Higgs phenomenology, it should satisfy two important ingredients for the no-lose theorem to 
discover a neutral Higgs boson at the Linear Collider (LC), the future e^e~ linear collider. 
One ingredient is that the lightest neutral Higgs boson should have a mass controlled by the 
electroweak scale, and the other is that at least one of the couplings of the Higgs bosons to 
a pair of Z boson is as large as that of the SM. 

In Ref. 12, Espinosa and Quiros investigated the Higgs sector of the (M+n)SSM, under 
the assumption of CP conservation. They calculated the upper bound on the lightest scalar 
Higgs boson mass in the (M+n)SSM by using the renormalization group approach. Accord- 
ing to their calculations, the lightest scalar Higgs boson mass of the (M+n)SSM is smaller 
than 145 GeV, if the gauge couplings in the theory are assumed to remain perturbative up 
to a very large energy, which they take as the Grand Unified Theory scale A = 10^^ GeV. 
The upper bound on the lightest scalar Higgs boson mass strongly depends on the scale A. 
In fact, they found that the upper bound on the lightest scalar Higgs boson mass increase 
from 145 GeV to 400 GeV, as A changed from 10^^ GeV down to 10^ GeV. 

We are motivated by the previous works of Ref. 12, where CP symmetry is conserved 
in the Higgs sector of the (M+n)SSM. In this article, we are interested in the possibility 
of CP violation in the Higgs sector of the (M+n)SSM. We would like to study whether the 
(M+n)SSM may accommodate complex phases in its Higgs sector for CP violation, while 
examining if the two above-mentioned phenomenological ingredients are, indeed, satisfied by 
it. An absolute upper bound on the lightest neutral Higgs boson mass will be obtained, and 
the possibility of discovering the lightest neutral Higgs boson at the LC will be discussed. 

The purpose of this article is not only to satisfy our theoretical curiosity about the 
(M+n)SSM but also to establish the phenomenological plausibility of the (M+n)SSM within 



the context of CP violation. In this article, we study CP violation in the Higgs sector of the 
(M+n)SSM. By taking the vacuum expectation values (VEVs) of the neutral Higgs fields 
and the relevant parameters in the tree-level Higgs potential as real, we assume that at the 
tree level, the Higgs potential of the (M+n)SSM has no complex phase. CP symmetry may 
then be violated at the one-loop level in an explicit way as radiative corrections are taken 
into account. For radiative corrections, we consider only the contributions from the third 
generation of quarks and squarks. 

We obtain the upper bound on the lightest neutral Higgs boson mass in the (M+n)SSM 
for given n by allowing the relevant parameters to vary within reasonable ranges; The mass 
depends evidently on n: The mass is found to increase as n grows to 10, but does not 
exceed 220 GeV for larger n. This maximum value of the upper bound is constrained by the 
condition that the squared masses of squarks be positive. Thus, production of the lightest 
neutral Higgs boson in the (M+n)SSM is kinematically allowed at the LC We, furthermore, 
investigate prospects for discovering the neutral Higgs bosons of the (M-l-n)SSM at the 
LC. The minimum cross section for producing any one of the neutral Higgs bosons of the 
(M+n)SSM in e^e~ collisions via the Higgsstrahlung process is found to decrease as n 
increases. 

This article is organized as follows: In the next section, we describe the Higgs potential 
of the (M+n)SSM and calculate the masses of the neutral Higgs bosons in the explicit CP 
violation scenario at the one-loop level. We then calculate the production cross section of the 
neutral Higgs bosons via the Higgsstrahlung process at the LC. Conclusions are presented 
in the last section. 

II. THE HIGGS SECTOR 

The Higgs sector of the (M+n)SSM contains the two Higgs doublet superfields TYf = 
{7ii,7i~) and Tif^ = (Ti^,?^^) of the MSSM and an additional n neutral Higgs singlet 
superfields A/^ (/ = 1 to n). In terms of these Higgs superfields, the most general form of the 
superpotential of the (M+n)SSM may be written as 

n n J 



where Aj and kiji {i,j,l = 1 to n) are dimensionless coupling constants, and e is the usual 
antisymmetric 2x2 matrix with ei2 = —£21 = 1. Here, for simplicity, the Yukawa couplings 
for the third generation are taken into account: The SU(2) doublet superfield Q^ = {ti, b^) 
consists of the left-handed quark superfields, the SU(2) singlet superfields t^ and 6^ are the 
charge conjugate of the right-handed quark superfields of the third generation, and hf and 
hb are the top and the bottom Yukawa coupling constants, respectively. 

At the tree level, the relevant Higgs potential of the (M+n)SSM consists of D terms, F 
terms, and soft terms as 

V^ = Vd + Vf + Vs , (1) 

with 
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J2 -fA,^^,N,N,N, + H.c. , (2) 

\i,j,i=i J 

where gi and g2 are the U(l) and the SU(2) gauge coupling constants, respectively, a = 

(or\ 0-^0-3) are the Pauh matrices, and Hf = {H^^H'), H^ = {H+,H^), and A^^ (/ = 1 to 

n) are, respectively, Higgs doublet fields and Higgs singlet fields. Note that the soft SUSY 

breaking masses rriHi, fnH2, and rriNi (i = 1 to n) are introduced in Vs, where A\- and A^. ,, 

{i,j, I = 1 to n) are the trilinear soft SUSY breaking parameters with a dimension of mass. 

We assume that the parameters in V"*^ are real. Thus, the CP symmetry is assumed to be 

conserved at the tree level. 

The total number of degrees of freedom in the Higgs sector of the (M+n)SSM from two 

Higgs doublets and n Higgs singlets is 2n+8. After the spontaneous breaking of electroweak 

symmetry, three of them are gauged away, becoming a neutral Goldstone boson and a pair 

of charged ones. The remaining 2n+5 degrees of freedom correspond to 2n+3 neutral Higgs 

bosons and a pair of charged Higgs bosons. In terms of the physical Higgs fields, after the 



spontaneous breaking of electroweak symmetry, the Higgs doublets and the n Higgs singlets 
may be expressed as 

, vi + hi+ i sin I3h^ 
Hi 

sin/3C+* 

COS/3C+ \ (3) 

-n2 = I 

V2 + /i2 + i COS I3h^ 
Nj = [xj + /i(2j+2) + ih(2j+3) ) , j = 1 to n , 

where hi {I = 1 to 2n+3) are the neutral Higgs fields, C"*" is the charged Higgs field, and 
Vi, V2, Xj (j = 1 to n) are the VEVs, respectively, of Hi, H2, Nj (j = 1 to n). We assume 
that the VEVs are real. With the above form for the Higgs doublets, one may take the 
constraints m'^ = g'^v'^ /2 and m| = [gl + gl^v^ jl for the electroweak gauge bosons with 



V = yv\ + V2 = 175 GeV. Note that tan/5 is given by the ratio of the VEVs of the two 
Higgs doublets, fi/f2- 

Now, we consider higher-order corrections. The Higgs sector of supersymmetric models 
at the tree level is well observed to be significantly affected by radiative corrections. The 
total Higgs potential at the one-loop level may be written as 

where V^ represents the one-loop corrections. The effective potential method allows us to 
obtain V'^ as [13] 

where A1| (/ = 1 to 4) and A4^ (g = b,t) are, respectively, the field-dependent squark and 
quark masses for the third generation, and A is the renormalization scale in the modified 
minimal subtraction (MS) scheme. We consider here the radiative corrections due to the 
loops of quarks and squarks of only the third generation. 

Note that the squark terms are positive whereas the quark terms are negative in V^. 
Because of the opposite sign between the squark and the quark terms, they would exactly 
cancel out if the masses of the quarks and the squarks were equal. Then, there would be no 
contribution to V'^. However, this is not the case because no squarks with masses comparable 
to quarks have been observed. Thus, we assume that the quarks and the squarks of the third 
generation have different masses, so the radiative corrections due to them should persist. 
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The masses of the stop and sbottom quarks are obtained as 

2 2 2|-'-/2i 2\-^ Pv^ 

"^i' ^t = ^b + T^i^Q + ^b) T \P^y , (4) 



2 

where mg, m^, and vhb are the soft SUSY-breaking masses, and vrit = htV2 and m^, = hi,vi 
are, respectively, the top quark mass and the bottom quark mass after the electroweak 
symmetry breakdown. The mixing between the two stop quarks is represented by X^, and 
the mixing between the two sbottom quarks by X-^. The two stop quarks or the two sbottom 
quarks will be degenerate in mass if Xi or Xi vanishes. We assume that they are not 
degenerate. They are, respectively, given as 

-. n n 

Xi = -{rriQ-m'^)^ + mf{A'f + {^X'jx^j)cot'^ [3 + 2 At{^ XjXj) cot [3 COS (l)t} , 

i=i i=i 

_. n n 

Xi = -{ml-mlf + ml{Al + C^ \^jx'^j)tiin'^ P + 2AbC^\jXj) tan p COS (pb} , (5) 

i=i i=i 

where At and A^ are the trilinear soft SUSY breaking parameters with a dimension of mass. 

We assume that At and Ab are complex, and we denote their phases by (pt and (pb, 
respectively. If At and Ab are real, (pt and (pb will be absent in the above expressions. Their 
presence implies CP violation at the one-loop level. CP symmetry is assumed to be conserved 
at the tree level, as the VEVs of the neutral Higgs fields and the relevant parameters in V'^ 
are taken as real. Now, CP symmetry is allowed to be broken explicitly at the one-loop level 
in the Higgs sector of the (M+n)SSM by complex At and Ab. 

If CP symmetry is conserved at the tree level and is broken explicitly at the one-loop 
level in the (M+n)SSM, there are non-trivial CP-odd tadpole minimum conditions for the 
neutral imaginary parts of Hi, H2, and Nj {j = 1 to n), namely, h^ and /i(2j+3) (j = 1 to 
n). Thus, the number of CP-odd tadpole minimum conditions in the (M-l-n)SSM is n -|- 1. 
They read 



3^f A(EJ=i ^jXj) sin0t , 2 2N ^mlAbiY^'-^i XjXj) sin (pb 
levrVsin^/? ■ly'^h^'^ij^ IGttVcos^/^ 



for h^ and 

SmfAtXjXj sin (pt ,, 2 2^ , ^mlAbXjXj sin (pb ,, 2 2n 
IQn^vsin^f] '' ^ '' '''' IQ-k^v cos^ {3 '^^ ^1' ^^^ 



for h2j+3 [j = 1 to n), where the scale-dependent function / is given by 

/("- '"?) = ji;^rh^^ {'"^ log § ' < log ^} + 1 ■ 

As can be seen from the structure of the above equations, the n+1 CP-odd tadpole minimum 
conditions may be reduced into just one equation, namely. 



m^Atsmcpt,, 2 2n , mlAbSmcpb 2 2^ 
sin2/5 ■'^ *i' *2^ cos2/3 •'^ ^1' ^^' 



sm^(3 ^^"^^^^ '^t,)^ ^^^2p /l^b'^. 
Further, the soft SUSY-breaking masses m'jj^, m'jj^, and mj^, {j = 1 to n) are eliminated by 
using the n+2 minimum conditions. 

The squared masses of the 2n+3 neutral Higgs bosons, ml (/ = 1 to 2n+3), are given by 
the eigenvalues of a mass matrix. These eigenvalues are obtained by the second derivatives 
of the Higgs potential with respect to the neutral Higgs fields hi {I = 1 to 2n+3), where 
the minimum conditions are employed. The symmetric (2ri+3)x (2n+3) mass matrix for 
the neutral Higgs bosons, M^^-, may be expressed in the basis of (hi, /i2, ■ ■ ■, /i(2n+3))- It is 
understood that, after diagonalizing the mass matrix, we sort the eigenvalues in increasing 
order by rearranging the neutral Higgs fields such that the (l,l)-element of the diagonalized 
mass matrix is Tn\^, the (2,2)-one is Tn\^, and so on. Thus, we implicitly assume hereafter 
that hi is not the neutral real component of the Hi Higgs field, but the lightest neutral 
Higgs boson, /i2 is not that of H2 but the second lightest neutral Higgs boson, and so on. 

Let us concentrate on the upper- left 2x2 sub-matrix of M?-. It is mathematically known 
that the upper bound on m\^ can be obtained from the upper-left 2x2 sub- matrix of M?-, 
regardless of the remaining matrix elements. The upper bound on m\ is then obtained as 
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+ f?^logHV^l (6) 

by using the two inequalities 

n n 

3=1 i=i 

n n n 

where 

n n 

^h = ^t COS 0t + ( J]^ Aj) ( J^ Xj) cot /9 , 

n n 

At2 = At + C^\j){^Xj) cot (3 cos (j)t , 
j=i i=i 

A^^ = Ab + (^\j)(^ Xj) tan P cos (pb , 

n n 

Al^ = AbCos0b + (^A,)(^x,)tan/5, (7) 

with the scale-independent function g given by 
(?K, ml) = -\ 1 log ^ + 2 . 



It can be observed that m^^ ^^^^^ depends on \j ( j = 1 to n) . As A/ increases, m^^ ^^^^^^ 
increases. It is worthwhile noticing that, if we replace YTj=i'^j by A and Yll=i^j by a;, 
the above expression for m\ ^^^^ looks exactly the same as the corresponding formula for 
the upper bound on the lightest Higgs boson mass in the (M+1)SSM, where A and x are, 
respectively, the dimensionless coupling constant and the VEV of the Higgs singlet in the 
(M+1)SSM. 

Now, in order to obtain the numerical value for m\^ ^^^, one has to take the values of \j 
and Xj {j = 1 to n), as well as other relevant parameter values. For simplicity, we set At = Ab 
and rriQ = tjit = mB- We choose a region of parameter space defined by 2 < tan/? < 40, 
< Aj < 0.7, < Xj,mQ,At < 1000 GeV (j = 1 to n), and tt/IOO < 0^,06 < 997r/100. 
We examine randomly 10^ points in the defined region of the parameter space. For each 
point, we check if rriq^ > rrit = 175 GeV and if rriq^ ^ rriq^ (g = t,6). The first condition 
is physically justified because it makes m\^ ^^^ more restrictive; and the second condition 



is required for numerical stability. If these two conditions are satisfied, we then calculate 
m\ ^^^. In this way, we obtain the upper bound on the mass of the lightest neutral Higgs 
boson in the (M+n)SSM for given n, which is an absolute upper bound in the sense that we 
examine the whole meaningful parameter space. 

Figure 1 shows m^ ^^^^ plotted against n, the number of Higgs singlets in the (M+n)SSM. 
As n increases, ml ^„^ increases from 148 GeV for n = 1. However, notice that it does 
not grow indefinitely. We find that it reaches 220 GeV for n= 15 and then very gradually 
decreases for n > 15. Thus, the universal upper bound on rrih-^ in the (M+n)SSM at the 
one-loop level for arbitrary number of Higgs singlets might be said to be 220 GeV. 

The upper bound on mh^ is actually set by the condition that the squared masses of the 
scalar quarks should be positive. From the expressions for the squark masses shown in Eq. 
4, it is clear that a/^ (g = t or b) might be restricted from above in order for the squared 
masses of the scalar quarks to be positive. Since the maximum value of itlq (= itlt = itlb) 
is on a SUSY-breaking scale (~ 1000 GeV), ^/X^ (g = t or b) must be smaller than (1000 
GeV)^. Otherwise, the electroweak symmetry breaking will not occur. Moreover, we have 
assumed that the lighter scalar quark masses are larger than the top quark mass. 

We note that in the (M+n)SSM, one can impose two additional constraints on the pa- 
rameter space. For any SUSY models in general, the heavier scalar quark masses must be 
smaller than the SUSY-breaking scale (~ 1000 GeV). That is, an absolute difference between 
the mass of any ordinary particle and the mass of its corresponding superpartner must be 
controlled by the SUSY-breaking scale. In particular, in SUSY models with Higgs singlets, 
the value of (X]?=i ^i) (Sj=i ^j) i^ limited by the electroweak scale rriEw ~ 300 GeV. If these 
two conditions are further imposed on the parameter space, m^ ^^^ should be lowered from 
220 GeV. In our calculations, these two constraints are not imposed. Thus, our result is a 
more conservative value. 

For the other neutral Higgs bosons, one can obtain an upper bound on ml (Z = 2 to 
2n+3) as functions of m.\^^^^ and m^^. Let O be the (2n+3)x(2n+3) orthogonal matrix 
that diagonalizes the neutral Higgs boson mass matrix. We now introduce 

Ri = Oil cos (3 + Oi2 sin P 

(/ = 1 to 2n + 3). It is clear that < i?^ < 1 and that a kind of orthogonality condition, 
namely, Yli=i -^f ~ l) i^ satisfied. The Ri represent collectively the whole parameter 
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space of the (M+n)SSM. Thus, it is more convenient to consider them than examining the 
parameter space itself because a particular set of Ri {I = 1 to 2n + 3) corresponds to a 
particular set of parameter values of the (M+n)SSM. In terms of Ri {I = 1 to 2n + 3), the 
upper bounds on the squared masses for the rest of neutral Higgs bosons are known to be 
given by 

Tfli, 



■•- {2^j=l^j) 

for / = 2 to 2n+3. 

In the (M+n)SSM, we expect that ZZhi (/ = 1 to 2n + 3), the coupling of a neutral Higgs 
boson to a pair of neutral gauge bosons, to decrease if n increases because the dimension of 
the neutral Higgs boson mass matrix increases by two for each additional Higgs singlet. Thus, 
for example, the production cross section of a neutral Higgs boson via the Higgsstrahlung 
process e^e~ -^ Zhi in e"'"e~ collisions decreases if n increases. In fact, the production 
cross section via the Higgsstrahlung process gets smaller if more Higgs singlets are added to 
the model (a large number of singlets is predicted in some string models) or if there is CP 
violation with a mixing of scalar and pseudoscalar Higgs bosons [14] . 

Now, let us consider the possibility of discovering any one of the neutral Higgs bosons 
in the (M+n)SSM at the LC We calculate the production cross sections of neutral Higgs 
bosons via the Higgsstrahlung process in e^e~ collisions. For the center of mass energy of 
the future LC, we take v^ = 500 GeV, 1000 GeV, and 2000 GeV. 

Whether or not any one of hi {I = 1 to 2n + 3) might be produced via the Higgsstrahlung 
process with real Z is determined by the threshold energy of Et = rnz + hi. For the LC with 
^/s = 500, 1000, and 2000 GeV, hi production is kinematically allowed because mh^ < 220 
GeV. If the production cross section of hi via the Higgsstrahlung process is small, one should 
then examine whether any of the other hi {I = 2 to 2n+3) has a production cross section 
large enough to be detected at the LC. 

Let o"sm(^ti) be the production cross section for the SM Higgs boson of mass m via 
the Higgsstrahlung process at a given center of mass energy. Then, in the (M+n)SSM, 
the production cross section a/ (/ = 1 to 2n+3) for the l-th neutral Higgs boson via the 
Higgsstrahlung process at a given center of mass energy is given by 

M^hi) = (rsM{mhi)Rf . 
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In particular, the production cross section for the heaviest neutral Higgs boson may be 
written as 

(2n+2) 



ff(2n+3)("^h(2„+3)) = ^Sm("^/.(2„+3)) 1 - J2 ^^ 



1=1 



Since the production cross section via the Higgsstrahlung process decreases as the neutral 
Higgs boson mass increases, we have for any of them 

a{m) > o-(mmax) = o-min • 

Hence, for given ruhi and Ri {I = \ to 2n+3), ai{min) satisfy 

O-i = 0-i(i?i,m/jJ > 0-l(m,j^_max) = CTl.min , 
(72 = a2{Rl,R2,mh.,) > 0-2(mft2,max) = 0^2 



min ) 



and so on. The minimum values, a^^min (/ = 1 to 2n+3), depend on m/i^^max and Ri {I = 1 
to 2n+3) such that 

0'2,min = ^2(^^/12, max) = 0'2,mm(-Rl) ""^fei.maxj j 

0"3,mm = <73("^/i3,max) = 0'3,mm(-Rl5 -^2, ""^/ii.max) 5 (8) 

and so on, because m| ^^^^ depend on /?/(/ = ! to 2nH-3). Examining these 2n+3 values 
o"«,min (/ = 1 to 2n+3), we may determine the possibility of detecting one of the 2n+3 Higgs 
bosons in e+e^ collisions. 

For given n, TTi/j^^max is fixed between 148 and 220 GeV. Then, for each n and for given 
R\ between and 1, we randomly generate i?f (/ = 2 to 2n+3) between and 1 with the 
constraint Yli=i -Rf = 1 in order to calculate o";,min (/ = 1 to 2n+3). Actually, this job 
is equivalent to exploring the whole parameter space of the (M+n)SSM. Among the 2n+3 
o"«,min (/ = 1 to 2n+3) for given n, let us select the largest one and denote it as 



a 



^rnax/ d2 



{Rl) . (9) 



This value is the minimum production cross section, above which at least one of the 2n+3 
neutral Higgs bosons is produced via the Higgsstrahlung process in e^e~ collisions. If we 
select, on the other hand, the smallest of them, it is the minimum production cross section. 
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above which all of the 2n+3 neutral Higgs bosons are produced via the Higgsstrahlung 
process in e^e~ collisions. 

We first set n = 1, and check the consistency of our analysis with previous studies in the 
(M+1)SSM. In the (M+1)SSM, there is only one Higgs singlet and five neutral Higgs bosons. 
We plot in Fig. 2 o-°'^^(i??) in the (M+I)SSM with exphcit CP violation as a function of i??, 
for ^/s = 500, 1000, and 2000 GeV. Our results are consistent with previous studies. For 
< i?i < 1, that is, for the entire parameter space of the (M+1)SSM, the a™^^{Rl) seems 
to be well-behaved. 

From Fig. 2, we find that the o-'"^^(i??) in (M+1)SSM at v^s = 500 GeV is no smaller 
than 10 fb. This implies that at least one of the five neutral Higgs bosons in the (M+1)SSM 
might be produced in e+e^ collisions for a certain set of parameter values of the (M+1)SSM, 
with its production cross section being larger than 10 fb. If the cm. energy of the e'^e~ 
system increases, we find that the minimum value of cr^'^^{Rf) in (M+1)SSM decreases: At 
y/s = 1000 GeV, a°'^^(i?2) is no smaller than 2.5 fb, and at v^ = 2000 GeV, a'"^^(i?2) is 
no smaller than 0.6 fb. Therefore, we may expect the integrated luminosity for 50 events 
of detecting one of the five neutral Higgs bosons in the (M+1)SSM to be about 10 fb~^ for 
the LC with ^ = 500 GeV, about 40 fb^^ for the LC with ^ = 1000 GeV, and about 
166 fb^^ for the LC with ^/s = 2000 GeV if we assume that the efficiency of the machine is 
about 50%. 

Now, let us consider the (M+n)SSM for n > 1. One may introduce a dimensionless 
parameter p in order to investigate the size of the CP mixing between the scalar and the 
pseudoscalar Higgs bosons. This parameter is defined in the (M+n)SSM by 



p = (2n + 3) 



2n+3/ d2 d2 



m p2 p2 



(10) 



The range of p is from to 1 because the Ri {I = 1 to 2n+3) satisfy the condition of 
'^1=1 -Rf = 1. If p = 0, there is no explicit CP violation in the Higgs sector of the 
(M+n)SSM. On the other hand, if p = 1, CP symmetry is maximally violated by the mixing 
between the scalar and the pseudoscalar Higgs bosons. The maximal CP violation that leads 
to p = 1 takes place when Rf = R2 = R^ = ■ ■ ■ = l/(2n + 3). In the case of maximal CP 
violation, the neutral Higgs boson search is extremely difficult. Actually, this can be seen 
in Fig. 2 for (M+1)SSM, where a°'^^(i?2) becomes smallest at Rj ~ 1/5. In the (M+1)SSM 
where we have five neutral Higgs bosons and R^ = R^ = ■ ■ ■ = R^, the maximal CP violation 
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will take place when Rf = R2 = ■ ■ ■ = R^ = 1/5. We deduce from Fig. 2 that the minimuni 
of a'^^'^{Rl) occurs at Rf ~ 1/5 in the (M+1)SSM when the rest parameters Rj {j = 2 to 
5) also become incidentally about 0.2 during random generation between and 1. 

For each n > 1, we repeatly determine the minimum of 0"™^^(i?^) for ^/s = 500, 1000, and 
2000 GeV. Our result is shown in Fig. 3, where the minimum of o'™'^^(i?f ) is plotted against 
n for ^/s = 500, 1000, and 2000 GeV. We see that the minimum of ^^^^(i?^) decreases 
monotonically as n increases: This is because the coupling coefficient of the Higgs bosons 
to the Z boson pair, at the ZZhj vertex, decreases as n increases. For example, for ^/s = 
2000 GeV, the minimum of (t'^^^{RI) is as small as about 0.08 fb for n = 15. 

In Fig. 4, we plot the integrated luminosity for 50 events of Higgs production against n 
in units of fb~^ for ^/s = 500, 1000, and 2000 GeV, where the efficiency of the e'^e~ colhder 
is assumed to be about 50%. For example, for n = 15, the integrated luminosity should be 
about 1084 fb~ for the future LC with y/s = 2000 GeV in order to detect 50 events of Higgs 
production via the Higgsstrahlung process. 

III. CONCLUSIONS 

We study the (M+n)SSM, a non- minimal supersymmetric standard model with n Higgs 
singlets and examine the possibility of CP violation in the model. We assume the CP 
symmetry is intact at the tree level, but is explicitly violated at the one-loop level in the 
Higgs sector. CP violation is induced by the mixing between scalar and pseudoscalar Higgs 
bosons due to the radiative corrections that are contributed by quarks and squarks of the 
third generation. In the presence of CP violation at the one-loop level, we calculate the upper 
bound on the mass of the lightest of the 2n+3 neutral Higgs bosons and their productions 
via the Higgsstrahlung process in e'^e" collisions at ^/s = 500, 1000, and 2000 GeV. 

The results of our study for n = 1 are consistent with previous studies for the (M+1)SSM. 
Our calculations yield an upper bound on the mass of the lightest neutral Higgs boson for 
n = 1 of about 148 GeV. The production cross section for at least one of the five neutral 
Higgs bosons of the (M+1)SSM via the Higgsstrahlung process in e~^e~ collisions at ^/s = 
500 GeV, according to our calculations, is no smaller than 10 fb. 

We also study the model for other values of n > 1. The upper bound on the mass of 
lightest neutral Higgs boson increases up to 220 GeV for n = 15. This value of 220 GeV 
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might be a universal upper bound on the mass of the hghtest neutral Higgs boson at the 
one-loop level in the (M+n)SSM for any number of Higgs singlets. This result suggests that 
the Higgs sector of the (M+n)SSM is regulated, showing no divergence or uncontrollable 
behavior for large n. 

The possibility of discovering at least one of the 2n+3 neutral Higgs bosons in e^e~ 
collisions via the Higgsstrahlung process was examined by calculating the minimum of the 
production cross section. We find that the minimum of the production cross section decreases 
monotonically for ^/s = 500 GeV, from about 10 fb for n = 1 to about 2.5 fb for n = 15. In 
terms of the luminosity, we require higher luminosity for the future LC if we try to discover 
at least one of 2n+3 neutral Higgs bosons in the (M+n)SSM via the Higgsstrahlung process 
for larger n. Our calculations suggest that an integrated luminosity of about 10 fb~^ for n 
= 1 to about 80 fb^^ for n =15 is required in order to detect 50 events of Higgs production 
at a future LC with ^/s = 500 GeV via the Higgsstrahlung process. We anticipate that a 
future LC with sufficient luminosity will yield experimental clues that may tell us how many 
Higgs singlets may be accommodated in the supersymmetric standard model. 
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FIG. 1: The upper bound on the mass of the hghtest neutral Higgs boson against the number of 
Higgs singlets, n. 
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FIG. 2: A plot as a function of Rf of the production cross section o"™^^(i?f ), as defined in the text, 
for at least one of the five neutral Higgs bosons in the (M+1)SSM via the Higgsstrahlung process 
in e+e~ collisions with ^/s = 500, 1000, and 2000 GeV. The would-be produced neutral Higgs 
boson is chosen such that the minimum of its production cross section is larger than the minimum 
of the production cross section of any other Higgs bosons in the model. Note that the curves have 
minima at R^ ~ 1/5. 



19 



ao(fb) 



10 


D 




















- 


D 


n 






(M+n)SSM 












D 












1 


_ A 

! ♦ 


A 


A 


A 


□ 
A 


D 

□ 

A 

A 


n 

A 


n 

A 


LC500 

° ° " = = „ 

LC1000 

^ ^ A A 

^ A A 


-1 
10 


- 






♦ 


♦ 




♦ 


♦ 


LC2000 


-2 
10 





















1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

n 



FIG. 3: A plot against n of the minimum of a^^'^{Rl) at i?f = l/(2n+3) for the future LC with 
Vs = 500, 1000, and 2000 GeV. 
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FIG. 4: A plot against n of the required integrated luminosity for 50 events of Higgs production 
via the Higgsstrahlung process in the future LC with ^/s = 500, 1000, and 2000 GeV, where 50% 
machine efficiency is assumed. 
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